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ABSTRACT 

We present the study of a quantum Controlled-Controlled-Not gate, implemented in 
a chain of three nuclear spins weakly Ising interacting between all of them, that is, 
taking into account first and second neighbor spin interactions. This implementation 
is done using a single resonant 7r-pulse on the initial state of the system (digital and 
superposition). The fidelity parameter is used to determine the behavior of the CCN 
quantum gate as a function of the ratio of the second neighbor interaction coupling 
constant to the first neighbor interaction coupling constant (J'/ J). We found that 
for J' /J > 0.02 we can have a well defined CCN quantum gate. 
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1. Introduction 



There is not doubt that the discovery of the polynomial time solution of the prime 
factorization problem (Shor's algorithm [1]) and the fastest data base searching 
(Grover's algorithm [2]) by a quantum computer, and the importance of the imple- 
mentations of these algorithms in cryptography analysis [3] have made the study 
of the physical realization of a quantum computer a priority for many researcher 
[4,5]. A quantum computer works with qubits instead of bit, where a qubit is the 
superposition of two quantum levels of the system which are defined as |0) and |1), 
and a tensorial product of these states makes up a register. Although, one could 
say that there are already quantum computers with registers of few qubits [6] , these 
are still far too low to make serious studies (which may require at least 100-qubits 
registers, that is, at least a 2 100 dimensional Hilbert space). In addition, Much more 
technology development (or very neat idea) is needed to have a real powerful quan- 
tum computer. One model of solid state quantum computer which looks promising 
due to advances in single spin measurements [7] is that of a chain of nuclear spins 
quantum computer [4,8], where the Ising interaction between first neighbor is con- 
sidered. This interaction allows us to implement ideally this type of computer up to 
1000 qubits or more [9]. First neighbor interaction between spins also allows to im- 
plement a Controlled-Not (CN) quantum gate using a single 7r-pulse [10]. However, 
the implementation of a Controlled-Controlled-Not (CCN) quantum gate, which is 
needed in Shor's factorization, Grover's searching and Full- Adder algorithms [11], 
is not that simple. In addition, CCN quantum gate is also important since it has 
universality characteristic [12]. The CN and CCN gates are defined classically by 
the following two tables [13] 
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where a (called control bit) and b (called target bit) for CN (a,b and c for CCN) 
represent the input information, and a' and b' (a', b', and c' for CNN) represent the 
output information. In a CN gate, "b"' (target) changes if an only if "a" (control) 
has the value 1. In a CCN gate, "c"' changes if an only if "a" and "b" has the value 

1. The ideals CN and CCN quantum gates, operating on an arbitrary L-register 
\iL-i, «a> ■■■)*&)■■■) *cj ■■■) *o) (where ij = 0, 1), would be defined as 

CNabllL-l, ...,i a ,...,i b ,...,i c ,...,io) = \lL-l, ... ,« , *a, ••• ,*c, ■•• ,*o) 

and 

CCN abc \i L -i, . . . ,i a , . . . ,i b , . . . ,i c , . . . ,i ) = \ i c ®(i a -ib), ■■■ ,io) , 

where the operation © means summation module 2. In this paper we show that 
considering second neighbor interaction between spins in the one-dimensional chain 
of spins quantum computer, and we show that it is possible to implement a CCN 
quantum gate with just a single 7r-pulse. Of course, a CN quantum gate continue 
being represented by a single 7r-pulse. 

2. Equation of Motion 

Consider a one-dimensional chain of three equally spaced nuclear-spins system (spin 
one half) making an angle cos (9 = l/Vs with respect the z-component of the mag- 
netic field (chosen in this way to kill the dipole-dipole interaction between spins) 
and having an rf-magnetic field in the transversal plane. The magnetic field is given 
by 

B = (bcosut, — bsinujt, B(z)) , (1) 

where b is the amplitude of the rf- field, B(z) is the amplitude of the z-component 
of the magnetic field, co is the angular frequency of the rf-field (its phase has been 
chosen as zero for simplicity). So, the Hamiltonian of the system is given by 

2 10 

H = - ]T ^ ■ B k - 2 Jh m+i - 2J'h E 4^+2 . ( 2 ) 
fc=o fc=o k=0 

where fi^ represents the magnetic moment of the kth-nucleus which is given in terms 
of the nuclear spin as fi^ = ^7(-^f, if, 1%), being 7 the proton gyromagnetic ratio. 
Bk represents the magnetic field at the location of the kth-spin. The second term 
at the right side of (2) represents the first neighbor spin interaction, and the third 
term represents the second neighbor spin interaction. J and J' are the coupling 
constants for these interactions. This Hamiltonian can be written in the following 
way 

H = H + W , (3a) 

where Hq and W are given by 

#0 = -n { E + 2 J{jm + if id + 2j'i z Q r 2 \ (36) 

U=o J 
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and 
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w =-—y: e**I+ + , 

fe=0 L 



-iuit 



(3c) 



Here, = r yB(zk) is the Larmore frequency of the kth-spin, 17 = 76 is the Rabi's 
frequency, and Ij^ 1 = 1% ± il% represents the ascend operator (+) or the descend 
operator (-). The Hamiltonian Hq is diagonal on the basis {(^l^o)} with ij = 0, 1 
(zero for the ground state and one for the exited state), 

H \i 2 iiio) = E i2ilio \i 2 iiio) ■ (4a) 
The eigenvalues E^i^ are given by 

Ei 2il i = -||(-ir 2 u; 2 +(-) ll Wl +(-iru;o+J[(-l) Jo+il + (-l) il+l2 ]+(-l) lo+l2 «/ / | . 

(46) 

The term (3c) of the Hamiltonian (3a) allows to have a single spin transitions on 
the above eigenstates by choosing the proper resonant frequency, as shown in Figure 
1. For example, if we are interested in having the transition 1 110) < — ► |1H) since 
this one represents the CCN quantum gate operation, according with the table 2, 
the resonant frequency would be 

to = ujq — J — J' . (5) 

Of course, we could also have considered a CCN quantum gate as defined by the 
transition | Oil) < — ► 1 1 1 1) , by defining the order of the elements differently. 



To solve the Schrodinger equation 



i*-Qt = , (6) 
let us propose a solution of the form 

*w = E<?*(*)|fc> , (6) 

k=0 

where we have used decimal notation for the eigenstates in (4a), Ho\k) = Ek\k). 
Substituting (6) in (5), multiplying for the bra (m\, and using the orthogonality 
relation (m\k) = 5 m k, we get the following equation for the coefficients 

7 

ihC m = E m C m + Y,Ck(m\W\k) m = 0,...,7. (7) 

fc=0 

Now, using the following transformation 

C m = D m e- iE ^ n , (8) 

the fast oscillation term E m C m of Eq. (7) is removed (this is equivalent to going 
to the rotating frame of reference), and the following equation is gotten for the 
coefficients D m 

I 7 

iL> m = t E W mk D k e^t , ( 9a ) 
n k=0 
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where W m k denotes the matrix elements (m\W\k), and uj m k are defined as 



Em - Ek 

h 



(96) 



Eq. (9a) represents a set of sixteen real coupling ordinary differential equations 
which can be solved numerically, and where W m k. are the elements of the matrix 



(W) = ~ 
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(9c) 



where z is defined as z = e , and z* is its complex conjugated. 



3. Numerical Simulations 

To solve numerically (9a), we shall use similar values for the parameters as reference 
(9). So, in units of 2ir x MHz, we set the following values 

uj q = 100 , wi = 200 , u 2 = 400 , J = 5 , ft = 0.1 (10) 

The coupling constant J' is chosen with at least one order of magnitude less that J 
since in the chain of spins one expects that second neighbor contribution to be at 
least one order of magnitude weaker than first neighbor contribution. We consider 
digital initial state and superposition initial states, 



*(0) = |110) (10a) 



and 



+ I^ |100> + 3! |101) + i> 10> + ^ |m> 



(10b) 



to cover all the important aspects of the CNN quantum gate. In all our simulations 
the total probability, J2 |Cfc(*)| 2 > is conserved equal to one within a precision of 10~ 6 
(note that \Ck\ = |-Dfc|)- Figure 2 shows the behavior of Re Dq, Im Dq, Re Dj and 
Im D7 during a 7r-pulse (t = r = tt/U) for the digital initial state and with J' = 0.1. 
One can see clearly the transition 1 110) — ► defining the CCN quantum gate 

up to a global phase (vr/2). Figure 3a shows the behavior of the probabilities \Ck\ for 
k = 0, . . . , 7 during a 7r-pulse for the superposition initial state and with J' = 0.1 
and for the superposition initial condition. Figure 3b shows the behavior of the 
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expected value of the z-component of the spin, (If.)-, k=0,l,2, during a 7r-pulse with 
the initial superposition state and J' = 0.1. These expected values are given by 



<tf> = ^E(-i) fc i<?fc(f)i 2 , 

fc=0 

(ID = \{\Co\ 2 + |Ci| 2 - |C 2 | 2 - \C 3 \ 2 + |C 4 | 2 + \C 5 \ 2 - |C 6 | 2 - |C 7 | 2 } , (116) 

(12) =\ii\c k \ 2 -izm 2 . (He) 

Z k=0 k=A 

As one could expect, only the spin related with the first qubit (having subindex zero) 
changes its value at the end of the 7r-pulse. The behavior of the expected values 
of the x and y components of the spin, (i|) and (i|f) for k = 0, 1, 2, is shown on 
Figure 4 for the superposition initial state and for f = 0.1. These expected values 
are given by 



and 

ITS) = 

2 



(IS) = Re (C^Co + C* 3 C 2 + C* 5 C 4 + C* 7 C 6 J , (I y ) = Im[. . . j , (12a) 
(if) = i?e f C 2 *C + C 3 *Ci + C 6 *C 4 + C$C 5 ) , (J*) = Im ( . . . ) , (12b) 



and 



(if) = iie (c 4 *C + C|d + C 6 *C 2 + C 7 *C 3 ) , (if) = Im (...) . (12c) 

The fast oscillations appearing on this behavior are due to the time depending phases 
in Eq. (8) which have these terms and which are not cancelled out. 

For the case J' = 0, the resonant frequency for the transition 1 110) < ► |1H) and the 

transition |010) < — ► 1 01 1) coincides and both transitions appears in the dynamics 
of the system (simultaneous CN and CCN operations). Thus, the question which 
arises is: what is the minimum value of J' which allows to have a well define CCN 
quantum gate ? To answer this question, Fig. 5 (a) and (b) shows the behavior 
of the CN quantum gate with the resonant frequency of the CCN quantum gate 
(uj = ujq — J — J') during a 7r-pulse for several J' values. In addition, (c) and (d) 
show the fidelity [14], F = (VP expected] 1 ® '( r )) > associated to CCN quantum gate as 
a function of J'/ J. As one can see from these figures, with J' values of even two 
orders of magnitude (J' / J = 0.02) we can already have a very well defined CCN 
quantum gate. Therefore, it is very likely that a single pulse CCN quantum gate 
can be implemented in a chain of nuclear spin quantum computer. 



4. Conclusions and comments 

We presented the study of a Controlled-Controlled-Not (CCN) quantum gate im- 
plemented in a chain of nuclear spins weakly Ising interacting and with the use of 
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a single 7r-pulse. With study this gate with digital and superposition initial states 
and found an expected global phase for its definition, 

CCN = |000)(000| + |001)<001| +|010)<010| +|011)<011| +|100)<100| 
+ |101)(101| + i|110)(lll| +i|lll)<110| . 

(13) 

Using the fidelity parameter, we have seen that it is very likely to have the imple- 
mentation of a CCN quantum gate in the chain of nuclear spin quantum computer 
since the coupling constant to second neighbor interaction can be even two order of 
magnitude lower than the first neighbor interaction, and still having a well defined 
CCN quantum gate. Since the eigenvalues in (4a) depends on J', Eq. (4b), the 
detuning factor, defined as A = (E p — E m )/% — uo for the transition \p) < — ► \m) 
(decimal notation), will depends also on J'. Therefore, the normal 27rk-method 
used in reference (8) will have to change correspondently to be able to suppress 
non-resonant transitions in multiple-pulses algorithms. 
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Figure Captions 

Fig. 1 Energy levels and resonant frequencies. 

Fig. 2 CCN quantum gate with digital initial condition and with J 1 = 0.1. (l):Re Dq, 
(2):im D 6 , (3):Re D 7 , (4):Jm D 7 . 

Fig. 3 For the superposition initial state and for J' = 0.1, (a) Probabilities \Ck(t)\ 2 
for k = 0, ... ,7. (b) Expected values (0): (Ig), (1): (If), and (2): (J|). 

Fig. 4 For the superposition initial state and for J' = 0.1, expected values of the 
transversal components of the spin. 

Fig. 5 For the superposition initial condition, (a) Probability | C2 1 2 , and (b) Prob- 
ability |C 3 | 2 for (1): J' = 0.0, (2): J' = 0.02, (3): J' = 0.04, (4): J' = 0.06, (5): 
J' = 0.08, and (6): J' = 0.1. (c) Real and Imaginary parts of the Fidelity, (c) 
Modulus of Fidelity. 
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